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SUMMARY 

The operation of the BGRID computer program is described for generating 
block-structured grids. Examples are provided to illustrate the code input 
and output. The application of a fully implicit AF (approximation 
factorization )-based computer code, called TWINGB (Transonic WING), for 
solving the 3D transonic full potential equation in conservation form on 
block-structured grids is also discussed. 
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I. INTRUOUCTION 


The basic feature of BGRID is that it generates grids with a 
block-structure (reference 1). An important aspect of BGRID is its 
adaptability to complex configurations with multiple components. Linear grid 
generation equations are used for simplicity. The grid produced is verified 
by computing the Jacobian for each grid cell. A good grid will have positive 
Jacobian for all grid cells to ensure that the transformation between the 
physical and the computational spaces is one-to-one. The block-structured 
grids generated by BGRID have been successfully used with finite volume 
solution algorithms using SLOR (successive line over-relaxation) in 2D 
(reference 2) and 3D (reference 1). 

This study adapts the TWINGB (reference 3) , an fully implicit AF-based 
computer code to the block-structured grids for the solution of the 3D 
transonic full potential equation in conservation form. This report describes 
how to use the modified TWINGB. 

The block-structurted grid generation methodology in the BGRID code is 
sumnarized in Chapter II. Input data description and output interpretation 
for BGRID are presented in Chapters III and IV. BGRID program information is 
given in Chapter V. Applications of the TWINGB code with block-structured 
grids are discussed in chapter VI. Operation of TWINGB is described in 
Chapter VII. 



II; BLOCK-STRUCTURED GRID 


The construction of a suitable grid system for complex 3D configurations, 
such as a wing/body/nacelle, is a necessary step for computing the 
corresponding transonic full potential flow. Two approaches have been 
available based on Thompson's (reference 4) surface-adapted coordinate 
concept. A limited approach maps the flow domain surrounding a 3D 
configuration into a single rectangular box. This approach has been 
successful for simple geometries (reference 5) but cannot be effectively 
applied to complex configurations with multi -components. A more general 
approach divides the computational domain into multiple rectangular blocks 
where the configuration itself is also represented by *a set of blocks whose 
structure follows the nature lines of the configuration (Figure 1). This 
results in a multi-block grid system which is adaptable to complex 
configurations with multiple components and can produce good grid quality near 
physical corners. A comparison of typical single block and multi-block grids 
is given in Figure 3. In General, a division of the flow field into multiple 
blocks allows the accommodation of slope discontinuities of the boundary 
surfaces and provide good boundary fitting behavior (Figure 4). 

A typical block-structured grid generation (Figure 5) process is described 
below: 

1. Define the overall block structure according to the natural lines of 
the configuration. 

2. Generate ID grids along the block perimeter (perimeter 
discretization). 

3. Generate 2D grids covering the block surfaces. 

4. Generate 3D grid covering the volume grids filling each block using 
the information obtained in Step 3 as a boundary condition for the 3D 
grid generation process. 

A simplified procedure is used to replace step 4 for wing grid 
generation. In step 3, 2Dv grids are generated for each wing span station. 
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Figure 2 , A Typical Block-Structured Grid. 
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Figure 4. Comparison of Grid Structure. 
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Figure 5. Block-Structured Grid Generation Process. 
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Connecting the cooresponding grid points between all neighboring span stations 
then produces a 3D grid. This code further assumes constant wing section so 
that the same 2D grid can be used for all span stations. 

Thompson’s concept of using a boundary-fitted curvilinear coordinate 
system requires that the curvilinear coordinates be solutions of a system of 
elliptic partial differential equations in the physical space subject to 
Dirichlet boundary conditions on all boundaries. One (curvilinear) coordinate 
is specified to be constant on each of the boundaries and a monotonic 
variation of the other coordinate along each boundary is specified. By 
interchanging the role of the Cartesian coordinates and the transformed 
coordinates, one arrives at a quasi-linear elliptic system for the Cartesian 
coordinates in the transformed space. 

The present approach simplifies the formulation by requiring the Cartesian 
coordinates x = (x, y) to be the solutions of three independent linear 
equations defined in the computational coordinates ( C > n): 

Ax__ + Bx + Dx_ + Ex =0 (1) 

CC nn 5 n 

where the coefficients A to E are constants or specified functions used for 
grid control. The coefficients A and B provide means to locally rescale the 
computational coordinates. Along the block boundaries, these two coefficients 
may be evaluated knowing the grid distributions. The coefficients D and E are 
related to the source terms in Thompson's nonlinear approach (reference 4). 
These two coefficients can be extracted effectively from the perimeter 
discretization of the block boundary (reference 6) by taking the limiting form 
of equation (1). For example, D may be extracted from the relation. 

Ax^^ + Dx^ =0 (2) 

The coefficients D and E are linearly interpolated on a block boundary surface 
using the block perimeter values. The grid control coefficients A to E are 
then linearly interpolated inside a block to match up with the boundary 
surface values. 
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III. INPUT DATA CARDS FOR BGRID 


This section describes the preparation of input data cards for the grid 
generation program BGRID. Assuming a user wants to analyze flow over a wing, 
he should provide a set of parameters for a block-structured grid. The user 
should also define the planform and section geometries of the wing. More 
information about the input variables will now be discussed. 

The case title card is read in with a 10A8 FORMAT. A set of data should 
be prepared for the NAMELIST GRIDIN if the user wants to override the program 
defaults. The variable heading cards are included for clarification purposes 
they may be replaced by blank cards if the user desires to do so. All the 
integer variables are read in with a 8110 FORAMT and all the floating point 
variables are read in with a 8F10.0 FORMAT. 
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Table 1 Data Deck for Blocked Grid Generation 


Columns 

Cards 

1-10 

11-20 

21-30 

31-40 

41-50 51-60 61-70 71-80 

Title of Case 


Data Cards for NAMELIST 

GRIDIN 

Title Card 1 

I BLOCK 

JBLOCK 

LINEAR 

NTHREE 

NPRINT 

Data Card 1 

4 

4 

1 

1 

1 

Title Card 2 

IBB 





Data Card 2 

1 

9 

31 

43 


Title Card 3 

JBB 





Data Card 3 

1 

9 

13 

21 


Title Card 4 

NBLOCK 





Data Card 4 

2 

2 

2 




2 

2 

-10 




2 

2 

2 



Title Card 5 

NCONTL 





Data Card 5 

4 

2 

2 




4 

1 

0 




3 

1 

1 



Title Card 6 

ILE 

ITE 

JLS 

JUS 


Data Card 6 

9 

33 

9 

13 


Tilte Card 7 

KMIN 

KTIP 

KMAX 



Data Card 7 

1 

11 

17 
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Table 1 (Continue) 


Tilte Card 8 

SWEEP 

ZTIP 

ZFAR 


Data Card 8 

30.0 

3.0 

9,0 


Tilte Card 9 

CHORD 

SPAN 



Data Card 9 

1.0 

0.0 




1.0 

0.1 




1.0 

0.2 




1.0 

0.3 




1.0 

0.4 




1.0 

0.5 




1.0 

0.6 




1.0 

0.7 




1.0 

0.8 




1.0 

0.9 




1.0 

1.0 



Title Card 10 

FSYM 

FNU 

FNL 

FMESH 

Data Card 10 

1.0 

72.0 

72.0 

1.0 

Title Card 11 

TRIAL 

SLOPT 

XSIN6 

YSING 

Data Card 11 

14.12 

0.0 

0.07 

0.0 

Title Card 12 

XB 

YB 



Data Card 12 

72 ^ 

cards 

for upper 

surface coordinates 

Title Card 13 

Not 

used 



Data Card 13 

Not 

used 



Title Card 14 

XX 




Data Card 14 

-10.0 

1.5 

1.5 

11.0 
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Table 1 (Continue) 


Title Card 15 

YY 




Data Card 15 

-10.0 

-0.3 0.3 

10.0 


Title Card 16 

DELI 

DELJ 



Data Card 16 

0.03 

0.015 
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Description of NAMELIST GRIDIN Parameters 


MAXIT Maximum number of iterations in the grid generation. Default value 
is 100. 

N0UT6 Interval of iterations for printout of transitory evolution of 
coordinates. Generally, the additional information will not be 
needed and the recommended value for N0UT6 is MAXIT+1. Default value 
is 101. 

IINCR Print the coordinates in the I-direction for every IINCR intervals. 
Default value is 4. 

JINCR Print the coordinates in the J-direction for every JINCR intervals. 
Default value is 4. 

LPLOT When generating the grid plot file, the program stores the 

coordinates for every LPLOT point in both the I- and J-directions. 
Default value is 1. 

OMEG The overrelaxation factor in the relaxation of the grid generation 
equations. Default value is 1.5 

ERR The convergence tolerance for the maximum corrections in the 
x-coordinates and the y-coordinates. Default value is 0.001. 

Data Card 1 


Column 

Variable 

Format 

Explanation 



1-10 

IBLOCK 

no 

Number of block boundaries 

•in 

I-direction. 

11-20 

JBLOCK 

no 

Number of block boundaries 

in 

J-direction. 

21-30 

LINEAR 

no 

Parameter to select grid generation equations 


1 use linear grid generation equations 
0 use nonlinear grid generation equations 
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31-40 NTHREE 110 Parameter to select 2D or 30 grids 

1 generate 3D wing grid 

0 generate 2D airfoil grid 

41-50 NPRNT no Parameter to control print out 

1 print convergence history 

2 print convergence history and block boundary 
grids 

3 print convergence history and block boundary 
and 2D grids 

Data Card 2 

Co 1 umn Variable Format Explanation 

1-11 IBB no Grid point indices of the block boundaries in 

. I-direction 

• 

71-80 

Data Card 3 

Column Variable Format Explanation 

1-11 JBB no Grid point indices of the block boundaries in 

. J-direction. 

71-80 


Example 
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JBB(4)=21 


JBB(3)=13 

JBB(2)»9 

JBB(l)al 


IBB(1)=1 IBB(2)=9 IBB(3)-33 IBB(4)=43 

This sketch illustrates how to use the IBB and JBB arrays to identify the 
I,J indices of the block boundaries in the global grid .systems. 

Data Card 4 

Co 1 umn Variable Format 

1-10 NBLOCK no 

71-80 
Example 


CM 

2 

2 

2 

1 

1 

o 

2 

2 

2 


Explanation 

NBLOCK denotes the block type I.O. number. The 
domain is filled with blocks and only some of 
them are filled with flow or wetted by the flow. 

1 block inside body (wetted by the flow) 

2 flow field block (filled by the flow) 

-10 block inside wake 
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This sketch illustrates typical block I.D. to be used for flow over an 
airfoil. 

Data Card 5 


Column Variable Format Explanation 






1-10 

NCONTL 

no 

Index to select the grid control parameters for 

• 



the interior grid points using block boundary 

• 



discretizations 

71-80 



0 linear variation between boundaries 




1 left and lower boundaries 



- 

2 left and upper boundaries 




3 right and lower boundaries 




4 right and upper boundaries 

Data Care 

1 6 



Column 

Variable 

Format 

Explanation 

1-10 

ILE 

no 

I-index for wing leading edge 

11-20 

ITE 

no 

I-index for wing trailing edge 

21-30 

JLS 

no 

J-index for wing lower surface 

31-40 

JUS 

no 

J-index for wing upper surface 

Data Card 

1 7 



Co 1 umn 

Variable 

Format 

Explanation 

1-10 

KMIN 

no 

K-index (spanwise) for wing root station 

11-20 

KTIP 

no 

K-index for wing tip station 
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21-30 KMAX 

Data Card 8 

no 

K-index for far field station in the spanwise 
direction 

Column 

Variable 

Forinat 

Explanation 

.1-10 

SWEEP 

FIO.O 

Wing sweep angle in degree 

11-20 

ZTIP 

FIO.O 

z-coordinate (spanwise) for wing tip station 

h-30 ZFAR 

Data Card 9 

FIO.O 

z-coordinate (spanwise) for far field station 

/ 

Column 

Variable 

Format 

Explanation 

1-10 

CHORD 

FIO.O 

Local chord length 

11-20 SPAN 

Data Card 10 

FIO.O 

Normalized spanwise coordinates 

Co 1 umn 

Variable 

Format 

Explanation 

1-10 

FSYM 

FIO.O 

Indicator to indicate whether or not the wing 
section is symmetrical 
0.0 nonsymmetry 
1.0 symmetry 

11-20 

FNU 

FIO.O 

Number of points to define the wing section 
upper surface 

21-30 

FNL 

FIO.O 

Number of points to define the wing section 
lower surface 

31-40 

FMESH 

FIO.O 

Use 1.0 
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Data Card 11 


Column 

Variable 

Format 

Explanation 

1-10 

TRAIL 

FIO.O 

Trailing edge closure angle 

11-20 

SLOPT 

FIO.O 

Trailing edge bisector slope 

21-30 

XSING 

FIO.O 

x-coordinate of the leading edge singular point, 
consistent with the wing-section defining data 
in data card 8 

31-40 

YSING 

FIO.O 

y-coordinate of the leading edge singular point, 


consistent with the wing-section defining data 
in data card 8 


Data Card 12 

Data card 12 contains as many cards as needed to define the wing upper 
surface section geometries. 


Column 

Variable 

Format Explanation 

1-10 

XB FIO.O 

x-coordinate of the wing section 

21-30 


defining points (upper surface) 

41-50 



61-70 




11-20 YB FIO.O y-coordinate of wing section 
31-40 defining points (upper surface) 

51-60 

71—80 

Data Card 13 

Data Card 13 should be included only if the corresponding wing section is 
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unsymmetrical (F5YM = 0.0 in data card 10). It should contain as many cards 
as needed to define the wing lower surface section geometries. 


Column Variable 

Format 

Explanation 

1-10 VAL 

FIO.O 

x-coordinates of the wing 

21-30 


section defining points 

41-50 

61-70 


(lower surface) 

11-20 DUM 

FIO.O 

y-coordinate of the wing 

31-40 


section defining points 

51-60 


(lower surface) 

71-80 


- 

Data Card 14 



Column Variable 

Format 

Explanation 

1-11 XX 

FIO.O 

Block boundary x-coordinates 

• 

• 

71-80 

Data Card 15 


at a far field boundary 

Column Variable 

Format 

Explanation 

1-11 YY 

FIO.O 

Block boundary y-coordi nates 

• 


at a far field boundary 


71-80 

Example 
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YY(4) 

YY(4) 

YY(4) 

YY(4) 



YY(4) 

YY(4) 

YY(4) 

YY(4) 


XX(1) • XX(2) XX(3) XX(4) 


This sketch illustrates how to use the XX and YY arrays to define the 
block boundary coordinates along far field boundaries. 

Data Card 16 


Column 

Variable 

Format 

Explanation 



1-10 

DELI 

FIO.O 

Minimum grid 

spacing in 

I-di recti on 

11-20 

DELJ 

FIO.O 

Minimum grid 

spacing in 

J-direction 
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IV. OUTPUT INTERPRETATION FOR BGRID 


The first page of output lists the parameters of NAMELIST GRIDIN. The 
second page prints the block structuring information. The third page prints , 
wing planform definition and the next three pages print the wing section 
defining geometries. The mapped coordinates are printed in the next two 
pages. Information about control parameters and grid convergence histories 
and coordinates are printed for each block in the following pages. An example 
for a block (16=4, JB=4) is demonstrated and grid information for other blocks 
is omitted because of its length. 

The grid is stored on TAPE4 for transonic potential flow calculations and 
a plot file of the grid is stored on TAPE? for graphical displays. Figure 6 
illustrates a. block-structured grid for a 3D wing. 
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V B6RID PROGRAM INFORMATION 


This chapter contains description of the major program modules of the 
computer program BGRID for blocked grid generation. B6RI0 consists of 15 
program modules: a main program (GRID) which contains the block-structured 

grid generation logic; several major subroutines which provide for block 
structuring, geometry specifications, grid initialization, grid relaxation, 
grid checking and grid output (SETUP, BOUND, MESH, MAPCO, STORE, DISC, 

INIGRID, CONTROL, RELGRID, SWEEP, GOUT); and three minor subroutines (INTPL, 
SPLIT, TRIMAT) which perform simpler functions, namely; interpolation, spline 
fitting, and solution of a system of tridiagonal equations. Each of these 
routines along with a brief description of its function is listed in table 2. 

A subroutine tree showing the relation of each routine to its called and 
calling routines is displyed in Figure 5. This information should be referred . 
to as needed during the following discussion. 

Description of MAIN Program GRID 

This is the driver program for wing block-structured grid generation using 
linear equations. The five subroutines, SETUP, BOUND, MESH, MAPCO and STORE 
are called by the main program sequentially. These five subroutines will be 
described in the following. 

Description of SUBROUTINE SETUP 

This subroutine reads in information for block-structuring. It also reads 
in geometric data for wing definition. The overall block-structuring should 
be set up by the user following the natural lines of the configuration. 

Description of SUBROUTINE BOUND 

This subroutine computes the block perimeters based on the information 
obtained from subroutine SETUP. The wing data is used to discretize wing 
surface using a parabolic transformation. SUBROUTINE DISC is called to 
discretize other block perimeter beside the wing surface. 

Description of SUBROUTINE MESH 


33 



This subroutine generates the computational grids in a block-by-block 
manner. For each block, SUBROUTINE INI6RID is called to provide an initial 
guess for the cartesian coordinates of the grid points. SUBROUTINE CONTROL is 
then called which computes the control parameters in the grid generation 
equations. Finally, the equations are solved through successive line 
overrelaxation (SLOR) by calling the SUBROUTINE RELGRID. 

Description of SUBROUTINE MAPCO 

This subroutine computes the Jacobian of each grid cell. MAPCO then 
checks to see if all 'the Jacobian are positive. Affirmative results indicates 
a useful grid with no cross-over. MAPCO will generate a warning message of 
BAD GRID-CROSS-OVER if any Jacobian becomes zero or negative. 

Description of SUBROUTINE STORE 

This subroutine saves the cartesian coordinates of the grid points in 
TAPE4 to be used by the TWINGB code. The subroutine also saves a plot file 
for the grid in TAPE7. 
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Table 2 Subroutine Reference Chart 


Routine 

name 

Calling 

routine 

Routines 

called 

Description 

BOUND 

GRID 

INTPL, SRLIF, 
EDGER ■ 

Calculates block perimeters 

CONTROL 

MESH 


Computes control parameters 

DISC 

BOUND 


Discretizes one block perimeter 

GOUT 

RELGRID 


Prints grids 

INIGRIO 

MESH 


Initiates grids 

INTPL 

BOUND 


Interpolation 

MAP CO 

GRID 


Checks the Jacobian of the grids 

MESH 

GRID 

INIGRID, CONTROL 
RELGRID 

Generates grids 

RELGRID 

MESH 

SWEEP, GOUT 

SLOR of grid equations 

SETUP 

GRID 


Reads block structure 

SPLIF 

BOUND 


Spline fitting 

STORE 

GRID 


Saves grid information 

SWEEP 

RELGRID 


Performs one sweep in the SLOR of 
grid equations 

TRIMAT 

SWEEP 


Tridiagonal solver 
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Figure 5. BGRID SUBROUTINE Tree 
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VI APPLICATION OF THE TWINGB CODE TO BLOCK-STRUCTURED GRIDS 


The computer program TWINGB (reference 3) uses a fast, fully implicit AF 
algorithm to solve the conservative full potential equation for transonic flow 
about arbitrary wings between walls. The basic TUINGB code uses single-block 
ring-grid. Current study adapts TWINGB to block-structured grids. 

In the ring-grid I » 1 and I = NI corresponding to the upper and lower 
surfaces of the vortex sheet, respectively. A periodic boundary condition is 
imposed to equate the values of the flow variables on the two sheets. In the 
block-strucutred grid 1=1 and I = NI corresponding to the upstream farfield 
and downstream farfield boundaries. Dirichlet boundary conditions are imposed 
along these farfield boundaries. 

Consideration of solutions accuracy requires that the mapping procedure 
between the physical and computational spaces be fully consistent. 
Specifically, the same approximation formulas should be used for the 
differencing of x,y ,2 and i in the computational space. The basic TWINGB uses 
second order formulas to difference i and uses fourth order formulas to 
difference x, y and z. This inconsistent procedure could provide satisfactory 
results for smoothly varying grids such as a simple ring-grid. But it will 
produce problems for general grids such as block-structured grids. The 
modified TWINGB uses the same second order formulas to difference x, y, z and 
i. 


Following the concept of isoparametric mapping, the difference formulas 
for x, y, z and i may be derived by differentiating t^e interpolators for x, 
y, z and-^. Same subdomain definition and interpolation formulas should be 
used for x, y, z and i. This fully consistent isoparametric mapping procedure 
willfurther improve the accuracy in the TWINGB flow solution algorithm but 
this procedure has not yet been implemented into the program. 

The complexity in the general block-structured grids calls for more 
program logics in the flow solver. Therefore, it is desirable to simplify the 
TWINGB code whenever possible. To maintain high computational efficiency 
including program vectorizability, it is necessary to remove the IF statement 
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from a OO-LOOP. To achieve this goal, calculations of flux quantities at 
different types of boundaries should be treated differently by a separate 
section of coding. This generates a large number of program instructions in 
the basic TWINGB code. Many of these instructions are only slightly different 
and repetition of nearly identical statements is error prone. In the modified 
TWINGB, these program instructions are packaged as statement functions. The 
use of nearly identical statements are eliminated by calling the required 
statement function. This also provides a better clarification of the 
program. Since the meaning of a statement function call can easily be 
associated with its name. An advantage in using the statement function over 
the conventional subroutine is that the statement function is an inline 
function. This means that at the compiling time, the effect of a statement 
function call is to replace the call statement by the corresponding program 
instructions which define the statement function. This would prevent the 
program from consuming CPU time in address transferring associated with a 
typical non-inline function call. Thus the modified TWINGB code maintains 
both computational efficiency and program simplification and clarification. 
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VII OPERATION OF TWINGB ON BLOCK-STRUCTURED GRIDS 


This chapter discusses the use of the modified TWINGB computer program on 
a block-structured grid including an input description for the program. Most 
of the parameter inputs to the modified TWINGB are contained in NAMELIST INP 
(defined in SUBROUTINE INPUT). These data include the flow parameters 
(freestream Mach number, angle of attack), and the aerodynamic and flow fiel 
convergence parameters. Default values for the NAMELIST parameters are 
intialized in INPUT. With the NAME LIST ’format, the user needs to change only 
the values that differ from the default values. The grid information is read 
in from TAPEIO. 

Description of NAMELIST INP Parameters 

ALPHAN Wing angle of attack in degree. Default value is 0. 

MACH Freestream Mach number. Default value is 0.72. 

ALPHAH Upper bound for the sequence of acceleration parameters in the 

AF scheme. Default value is 4.0. 

ALPHAL Lower bound for the sequence of acceleration parameters in the 

AF scheme. Default value is 0.40. 

BETA Temporal damping coefficient in the ETA direction. Default 

value is 0.0. 

BXI Temporal damping coefficient in the XI direction. Default value 

is 0.0. 

Maximum number of iterations allowed per case. No default value 
given. It must be supplied by users. 

Relaxation parameter in the iteration scheme. OMEGA should be 
bounded between I and 2 inclusively. Default value is 1.8. 


MAXIT 


OMEGA 
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CON 

M 

KKK 

BB 

S 

lOUT 

JOKIP 

INTOUT 


Parameter controlling upwind bias of density. CON should be 
bounded between 1 and 2 inclusively. Default value is 1.5. 

Number of elements in ALPHA sequence. M must be no less than 
2. Default value is 8. 

Starting value of ALPHA sequence. KKK may vary from 0 to M. 
Default value is 0.. 

Wing span calculated from the grid. 

Same as BB. 

Output option parameter. Should set a default value of 0 to 
print the density and Mach number and plot Cp. 

Skip parameter in space for output plots. 

Skip parameter in time for output plots. Printed after 
convergence only if set to 0. Default value is 0. 
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